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Abstract. The critical behaviour of d-dimensional semi-infinite systems with 
n-component order parameter <p is studied at an m-axial bulk Lifshitz point 
whose wave- vector instability is isotropic in an m-dimensional subspace of K''. 
Field-theoretic renormalization group methods are utilized to examine the special 
surface transition in the case where the m potential modulation axes, with 
< m < d— 1, are parallel to the surface. The resulting scaling laws for the surface 
critical indices are given. The surface critical exponent r;n^, the surface crossover 
exponent $ and related ones are determined to first order in e = 4+ ^ — d. Unlike 
the bulk critical exponents and the surface critical exponents of the ordinary 
transition, # is m dependent already at first order in e. The 0{e) term of rjn^ is 

found to vanish, which implies that the difference of fS'^ and the bulk exponent 
P is of order e^. 
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1. Introduction 

As is well known, systems in critical or near-critical states are sensitive to the 
presence of boundaries such as surfaces or walls: the long-distance behaviour of 
their local densities and correlation functions gets modified near boundaries. During 
the past decades impressive progress has been made in our understanding of such 
boundary critical phenomena (Binder 1983, Diehl 1986, Cardy 1987, PleimHng 2004). 
A prototype class of systems from whose study considerable insight into general 
aspects of boundary critical phenomena has emerged is provided by the semi-infinite 
n-vector cf>'^ models. At bulk criticality, these models exhibit a variety of physically 
distinct continuous surface transitions in sufficiently high space dimensions d, called 
ordinary, special, extraordinary, and normal. A simplifying feature they have is 
that the scale invariance they display at criticality is of isotropic nature: distances 
along arbitrary directions must be scaled with the same power of the length-rescaling 
factor i. Furthermore, the correlation lengths ^(e) defined through the exponential 

X On leave of absence from Institute of Solid State and Semiconductor Physics, 220072 Minsk, 
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decay ~ exp(— i?/^(e)) of the pair correlation function {(p{x) (p{x + Re)) along a 
given direction specified by a unit vector e all diverge ~ |T — TJ, d"" with the same 
critical index v, as the critical bulk temperature T = T^^h is approached via a phase 
with unbroken 0{n) symmetry. This leaves little room for dependence of boundary 
critical behaviour on the orientation of the surface with respect to crystal axes. We 
are aware of only one established scenario for such orientation dependence (Schmid 
1993, Diehl et al 1997, Leidl and Diehl 1998, Leidl et al 1998): bcc binary alloys 
and antiferromagnets in the presence of weak magnetic fields may map upon coarse 
graining onto semi-infinite (/)^ models with or without ordering boundary fields /ii, 
depending on whether their surface's orientation breaks the symmetry cf) -^ —(j) or 
not. If so, they may belong to basins of attraction of distinct fixed points, such as the 
normal and ordinary one, so that their boundary critical behaviour is different. 

In the present paper we shall be concerned with boundary critical behaviour at 
m-axial bulk Lifshitz points (LP) (Hornreich 1980, Selke 1992, Diehl 2002). An LP is 
a multi-critical point at which a disordered, a homogeneous ordered and a modulated 
ordered phase meet. A characteristic feature of the scale invariance that systems 
exhibit at such points is its anisotropic nature: the coordinates Xa, oi — l,...,m, 
along any of the m potential modulation axes scale as a nontrivial power {xp) of 
the d — m remaining orthogonal ones, xp, (3 — m + 1^ . . . ,d. This entails that the 
orientation of the surface with respect to the modulation axes matters. Two principal 
surface orientations, parallel and perpendicular, can be distinguished for which the 
surface normal is directed along an a or /? axis, respectively. The boundary critical 
behaviour that occurs at a given type of surface transition (ordinary, special, etc) for 
either parallel or else perpendicular surface orientation is expected to correspond to 
two distinct universality classes. Results obtained via mean-field theory (Gumbs 1986, 
Binder and Frisch 1991, Frisch et al 2000) and Monte Carlo simulations (Pleimling 
2002) for the ordinary and special transitions of ANNNI models with free surfaces of 
both types of orientation are in conformity with this expectation. 

These findings indicate that systems exhibiting anisotropic scale invariance have 
potentially richer boundary critical behaviour. Since systems of this kind are abundant 
in nature — anisotropic scale invariant states exist both in and out of equilibrium — a 
careful understanding of their boundary critical behaviour is certainly of interest. 

Focusing on the case of parallel surface orientation, we have recently introduced 
an appropriate semi-infinite extension of a standard n-vector (f)^ model describing 
universality classes of bulk critical behaviour at m-axial LP points (Diehl et al 
2003a, 2003b). We argued that this extension represents the surface universahty 
classes of the ordinary, special and extraordinary transitions for this kind of surface 
orientation (Diehl 1986, 1997). A renormalization group (RG) analysis for dimensions 
d = d* (m) — e below the upper critical dimension d* (m) = 4-1- to/2 lent support to this 
claim, giving fixed points we could identify as describing the corresponding surface 
universality classes. However, this identification was tentative: since of all, only the 
ordinary transition was analysed in some detail; second, this was done by relying on 
the asymptotic validity of the Dirichlet boundary condition, bypassing thereby the 
need to determine the precise location of the ordinary fixed point in the space of 
surface interaction constants. 

The purpose of the present paper is to complement our previous work (Diehl et al 
2003a, 2003b) by a detailed analysis of the special transition. The asymptotic (Robin- 
type) boundary condition that the theory satisfies at this transition involves the fixed- 
point value A!j_ of a dimensionless surface variable A, a number which depends on the 
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order of the loop expansion. A trick analogous to the one employed in the analysis 
of the ordinary transition by which the computation of X*^ can be avoided does not 
exist. This entails a qualitative difference between the analysis of the m > special 
transition described below and that of its to = counterpart at bulk critical points 
(Diehl and Dietrich 1981, 1983): in the latter case, the variable A is missing, and 
the fixed-point value of the sole remaining renormalized surface interaction constant 
(the surface enhancement c) is zero to all orders in perturbation theory. Unlike the 
TO. > case, projection onto the hyper-plane of critical renormalized surface variables 
is trivial. 

We shall determine the surface correlation index ry^^ and the surface crossover 
exponent $ to first order in e. Our RG analysis reveals that below the upper critical 
dimension, all other surface critical exponents of the special transition can be expressed 
in terms of these two independent surface critical exponents, along with bulk critical 
indices. While this is just as in the to = case of the special transition at bulk 
critical points (Diehl and Dietrich 1981, 1983, Diehl 1986), our result that the 0(e) 
contribution to m^ vanishes if to, > (though not for ?Ti = 0) is a remarkable difference. 
Another unusual feature of our results is that the crossover exponent $ starts to 
depend on m already at linear order in e. By contrast, the e-expansions of both 
the bulk critical exponents as well as the surface critical exponents of the ordinary 
transition display first ?Ti-dependent deviations from their m = analogues only at 
quadratic order in e. 

The remainder of this paper is organized as follows. In the next section, we define 
our model and briefly recapitulate the relevant background knowledge required in what 
follows. Section 3 deals with the RG analysis of the special transition. We derive the 
scaling forms of the multi-point cumulants involving order parameter fields (p on ^^^^d 
off the boundary, give some of the implied scaling laws, and present our e-expansion 
results. The latter are utilized to estimate the values of the surface critical exponents 
rh^ and /3^^, as well as the surface crossover exponent for the uniaxial scalar case in 
three dimensions, i.e. the ANNNI model. In section 4, RG improved perturbation 
theory is employed to investigate the surface susceptibility xii- Making an explicit 
one-loop calculation, we verify the predicted scaling behaviour, corroborating thereby 
our identification of the special transition. A short discussion and concluding remarks 
follow in the final section. Finally, there are two appendices explaining computational 
details. 

2. Model and background 

Following Diehl et al (2003a, 2003b), we consider a model defined by the Hamiltonian 
H= f Ch{x)dV+ f Ci (x) dA , (1) 

with the bulk density 

^a=l ^ /3=m+l a=l 

and the surface density 

o ? m 

A(a;)^^02^-^(a„0)\ (3) 

a=l 
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where the volume and surface integrals L, dV and J^ d^ extend over 
the half-space 5J = {a; = (r, z) | r e R''^^, < z < 00} and the boundary plane 
*B = {(r,0) I r e R'^"^} ('the surface'), respectively, and the notation da = d/dxa 
and 9/3 = d/dxp^ with 1 < a < m and m < (3 < d, is used. The order parameter is an 
n- vector 4> — {4>a, o, = 1, . . . , ri). 

Let us introduce the {N + M)-point cumulants involving N fields (/>a. (aJi) off the 
surface and M boundary fields 4>^.{rj) = (j)bj {rj,z = 0): 



N n r *^ n \ cum 



n-^a-c^o n<^^(^^) 



3= 



(4) 



where a; and r are shorthands for the sets of all points Xi and Vj off or on the 
boundary, respectively. From previous work (Diehl et al 2003a, 2003b) we know that 
the ultraviolet (uv) singularities of these functions can be absorbed for d < d* (m) by 
means of the following re-parametrizations: 

4> = Zj 0rcn , (7 = Z^a , u a^™-!^ F„i,^ ^ ^^ ZuU , 

T - TLP = H"^ Zr[T + Ar p^] , {p - php) CT^^^^ = f^ Zp p , (5) 

(/.'^ = (Z0Zi)i/VL , A (7-1/2 =x + Px{u,X,e), 

c~Csp ^ pZc[c + Ac{u, A, e) p] , (6) 



where 



r(l + e/2) r2(l - e/2) r(m/4) 



^"'■' (47r)(8+™-2^)/4r(2-e)r(m/2) ^''^ 

is a convenient normalization constant. Further, tlp and /5lp are the critical bare 
values of the LP; in our renormalization scheme based on dimensional regularization 
and the e-expansion, these fluctuation- induced shifts vanish. The same applies to Cgp, 
the critical value of the bare surface enhancement c at which the special transition 
occurs. 

The bulk renormalization factors Z^^r,p,u — Z^,r,p,u{u, e) are power series in u, the 
renormalized coupling constant, and Laurent series in e; since we employ the scheme 
of minimal subtraction of uv poles, their regular part at e = is exactly 1. The 
boundary renormalization factors Zi^c = Zi_c{u, A, e) have a similar form, except that 
the series coefficients depend on A. Finally, P\ and Ac are renormalization functions 
of the form 

oo oo oo 

P,{u,X,e) = J2 P^^-'^W^^' = E Y.P^'''''^^^'>^ ■ (8) 

i,j — l *,i — 1 k—0 

Two important features should be appreciated: first, although the 0{u) term of 
P{u, A, e) is proportional to A, its ©(u^) contribution is not, and hence does not vanish 
for A = 0. Thus P{u, A = 0, e) 7^ 0, i.e. a nonzero value of the bare interaction constant 
A gets generated by the 0^ interaction (Diehl et al 2003a). Second, the renormalization 
of f and c mixes the corresponding renormalized quantities r and c with p^ and p, 
respectively. The consequences (worked out in (Diehl et al (2003a) and below) should 
be no surprise: the role of the scaling fields r and c of the Gaussian theory is taken 
over by combinations of the form 

9t{t,p,u) = T + Cl2{u)p'^ , gc{c,p,u,X) =c + Cp{u,X)p . (9) 
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To become more specific, it is necessary to recall the RG equations of the 
renormahzed cumulants gS"^^ - ^-iN+M)/2^-M/2 ^(^nm) ^ They read 



p— 'U,(7",T,p,C,A 



N + M 



M 



■v<p- 



Vi 



GLr')=0. (10) 



The beta functions, defined by /3p = /i9^Lp, where 9p|o denotes a derivative at 
fixed bare interaction constants, can be expressed in terms of the exponent functions 
r?p,0,i = Mi9^|QlnZp,0,i and Pa as 

/3„(u, e) = -u[e + riu{u)] , 

/3cr(u,cr) = - crrj^{u) , 

f3r(u,T,p) = -T[2 + lJr{u)]-p^br{u) , (11) 

Pp{u,p) = - p[l + r]p{u)] , 

Pc{u, X,p,c)^ - c [1 + 77c(u, A)] - p bc{u, A) , 



and 



Px{u,\) = 



-(3u{u,e)duP\{u,\e) 



(12) 



l+dxPx{u,\,€) ' 

where the contributions involving b^. and be are produced by the terms proportional 
to At and Ac in equations (5) and (6), respectively. 

Explicit two- loop results for the bulk quantities Z(j,^a-,T,m Pu, and ry^.o-.r.p can 
be found in Shpot and Dichl (2001); see also Diehl and Shpot (2000). The bulk 
function A^- is given to order u in equation (17) of Dichl et al (2003a), where also the 
renormalization functions Zi, Zc, and P\ of the semi- infinite system were obtained 
to the same order in u for general values of A > 0; see its equations (21)-(23). The 
remaining surface counter-term Ac is computed in Appendix A; we obtain 



Ac{u,\,e) 



2 u 



l-n(A;m)- Ai'i(A;TO) V^r[(m -t- 2)/4)] 



2A 



2r(m/4) 



+ 0(^2). (13) 



Utilizing these results gives the beta function 

n + 2 



/3a(w, A) = -- 
and the exponent functions 

?7i(m,A) = -■ 
and 

r]c{u,\) = -- 



2 p(2,-l)/ 



6 



■ii(A;m)AM + 2M^P^"' '^(A)+0(w') 



■ ii(A; m) u + 0(u 



6 



[l + Xi[{X;m)]u + 0{u^) 



Here ii{X;m) and ii{X;m) are integrals defined through 
ii{X;m) = ii(A, e = 0;to) , 
i'i{X;m) = i'i{X, e = 0; m) , 
i'i{X, e; m) = dii{X, e; m)/dX , 

«i(A, e;m) 



1 + At 



where we have introduced the convenient notation (/(t))^ 



(14) 
(15) 
(16) 

(17) 
(18) 
(19) 

(20) 
for the normalized average 



{fit))e,r. 



K 



£,m Jo 



di/(t)t("-2)/2(^_^2)(2-2e-™)/4 



(21) 



^,(A;™).2.F,(l,f;|;A^)-l+V^^ii^^_Mi^ 1 - (l - A^)(-")/^l , (25) 



Boundary critical behaviour at Lifshitz points 6 

over t. The normalization factor is given by 

M,m ^ / d«("-2)/2 (^ _ ^2)(2-2e-m)/4 ^ 1 5[^/4^ (6 - 2e - m)/4] (22) 

Jo 
where B{a, b) is the Euler beta function. This choice of A/'e,m ensures that 

ii(0,e;m) = l. (23) 

The t-integral required for ii{X]m) converges for < to < 6. Since the 
norniahzation factor varies as l/A/o.m ^ (6 — to,) for m —^ 6, the integral's singularity 
^ (6 — m)^^ gets cancelled to produce the finite limiting value 

^l(A;6-)= lim zi(A;to) = 1— -. (24) 

m^6— i + A 

For general values m G (0,6), the integral ii(A;?Ti) can be expressed in terms 
of hyper-geometric functions 2-F1 and elementary ones. Term-wise integration of the 
integrand's Taylor series in A or evaluation via Mathematica§ leads to 

/^r[(m-2)/4] 
r(m/4) A 
which simplifies considerably if jti = 2 or to, = 6, giving 

zi(A;2) = 2A"Mn(l + A)-l (26) 

and 

TT 2 arccos(A) , , 

respectively. Explicit plots of the functions — Aji(A;to,) for m = 1,2,..., 6 are 
displayed in figure 1 of Diehl et al (2003a). 

Setting the relevant bulk variables t and p to zero and the coupling constant u 
to the nontrivial zero of /?„ for e > 0, namelyj 

u*^^^ + 0{e^), (28) 

3 n + 8 

we can analyse the RG flow and determine the fixed points in the cA plane. This 
yields the schematic flow picture shown in flgure 1 Diehl ct al (2003a). 
The fixed points 

^o*rd : (cSrd - «^, A = a;) , 

n*p: (c*p = 0, A = A;), (29) 

^ox : (ccx = -00, A = a;) , 



located at the positive value 

72 ] 

a;(to) = Ao(to) + ^— ^^^^-— ^i:^^^^:- - + o(6'), (30) 



72Pf-^)(Ao)e 
,+ ,».; - .ov'V ^ („ + 2)(n + 8) Ao ii(Ao; to) 



should describe the ordinary, special and extraordinary transitions, respectively. Here 
Ao = Ao(to) arc the real positive zeros of the functions «i(A;to,); for later use we list 
their values in table 1. Wc have included the limiting values for jti ^ 0. Of course, for 

§ Mathematica, version 4.1, a product of Wolfram Research. 

II The expansion of u* to 0{e'^) can be found in equations (60) of Shpot and Diehl (2001) or equation 
(42) of Diehl et al (2003a). It will not be needed in the following because our subsequent analysis is 
based on one-loop results for the surface RG functions fie and »yi,c- 
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A; = o(eO) 



y = o{e) 



Figure 1. Schematic picture of the RG flow in the c\ plane at p = r = and 
u = u* , showing the fixed points V*^^, 'Psp^ ^'^^ ^ox specified in equation (29). 



Table 1. Zeros Ao(m), derivatives i'^{\;m) and coefficients 1 + Aoij(Ao;m) 
appearing in the 0(e) terms of r]*; see equation (60). In the second column 
the exact limiting values for m — > are listed. The values for m = 1, . . . , 5 are 
approximate numerical numbers; those for m = 6 are exact. 



m 







1 


2 


3 


4 


5 


6 


Ao(m) 




oo 


6.20921 


2.51286 


1.70176 


1.34277 


1.13629 


1 


-i[[Xo{m);m] 







0.05505 


0.17138 


0.27332 


0.35949 


0.43400 


1 

2 


1 + \o{m) i[[Xo{m); 


to] 


1 


0.65816 


0.56934 


0.53488 


0.51729 


0.50685 


1 
2 



m = 0, the variable A becomes meaningless. However, considering to- as a continuous 
variable, we can ask what happens to the RG flow, the zero Ao(to.) and the critical 
exponents as to, ^ 0. Utilizing either the integral representation (17) of «i(A;to,) or 
the explicit form (25) one easily verifies that limm_>o«i(A; to) = 1 for arbitrary fixed 
A > 0. Hence the root Ao(to.) must approach infinity in this limit. To obtain its 
limiting behaviour, one can work out the asymptotic expansion of ii for A ^ oo at 
fixed TO, < 2. A straightforward analysis yields 



ii(A;TO) = 

A— *oo 

which iniphes the small-m form 



-l + A-/^J"t!!!^"^"('?.+0(A-^) 



r(TO/4]r[6 - to)/4] 



Ao(m) « 2(2-™)/™ e 

■m—*0 



(31) 



(32) 



An evident consequence of these two equations is that the product Ao(m) i'l [Ao(m); m] 
(which will be needed below) vanishes for to, ^ 0. 



^:7^ pW = /3p["W' ■ • ■] ' P(1) = P' P^u,a,T,p,c,X, (34) 
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3. RG analysis of the special transition 

In order to analyse the asymptotic behaviour described by the special fixed point 
P*p, we must determine the relevant bulk and scaling fields that can be constructed 
from T, p, and c. Following Diehl et al (2003a), we choose the coefficient C\{u) in 
equation (9) such that the running variable gr{i) into which g^- transforms under a 
scale change p ^ pi satisfies 

^Jg9Ai) = -[2 + Vr{u)]9r{e), 5r(l) -.9r, (33) 

where the asterisk on functions of u means their value at u = u* . 
Upon taking into account the flow equations 

d_ 

one easily derives the condition 

hr{u) = [rjr{u) - 2T^p{u) + Pu{u,€)du]ClAu) , (35) 

which is easily solved for u near u* to obtain 

^;^(") = -^ + 0(S-^.*). (36) 

'It '^'ip 

We fix the coefficient Cp of the scaling field Qc in a similar fashion. Requiring 
that 

i-^^9cie) = -[l + Vciu,\)]9c{'i), 5c(l)=5c, (37) 

we arrive at the condition 

6e(u,A)= [r^c{u,~X) - r^piu) + f3^iu,e) d^ + Pxiu,~X) d-^]C'p{u,~X) , (38) 

which in turn yields 

C,^(w,A) = ^^ + 0(s-u*,A~A;) (39) 

Vc 'Ip 

for (u, A) sufficiently close to the fixed-point values (u*,A^), where the asterisk on 
functions of u and A indicates their values for {u, A) = (u*, X*^), i.e. at the fixed point 
7?* 

' sp' 

Clearly, the mixing of r and c with p should vanish as m ^ 0. This is indeed the 
case as we see from the explicit 0(u) results 

, , , n + 2 mu „ , o, , , 

hr{u)^^ — + 0{u') (40) 

and 

"^ ' ' 3 |_ 4A 4r(m/4) J ^ V y ' ^ > 

the first of which was obtained in Dichl et al (2003a), while the second follows from 
equation (13). Note that the 0(u) terms of both 6t- and be vanish in the limit to ^ 
quite generally, not only at the fixed point V*. 

From equations (33) and (37) we can read off the RG eigenexponents of the scaling 
fields Qt and gc- Writing them as 1/v and ^ /v, respectively, we identify the standard 
correlation-length exponent 

z.= (2 + <)^i (42) 
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(denoted by 1^12 in Shpot and Dichl (2001), Diehl et al (2003a)) and the surface 
crossover exponent 

$ = z.(l + r;:). (43) 

Let us also introduce the bulk anisotropy exponent 

^=i(2 + r/:), (44) 

the familiar bulk critical exponents 

7,^77;, p^\u[d-2 + i^ + m{e-l)], (45) 

(denoted by 77/2 and /3;, respectively, in Shpot and Diehl (2001), Diehl et al (2003a)), 
the bulk crossover exponent 

^ = u{l+r^;) (46) 

ndices 

(47) 

^{9 - 1)] = /3 + \ui^l (48) 

of the special transition. 

It is now straightforward to solve the RG equation (10) via characteristics 
and derive the asymptotic scaling forms of the cuniulants G^^'^' near the special 
transition. Setting /i = cr = 1 for notational simplicity and neglecting corrections to 
scaling, we obtain 

Gien^^\ra,rp, z; p, gr, u, gc, A) 

The special choice iV = and M ~ \ shows that our identification (48) of (f^ is 
correct. On the other hand, our definition (47) of r/^J' may need explanation. Noting 

that the surface cumulant Gron (''jO) should have a well-defined limit as ^r, cr and 
gc approach zero, we can read off from equation (49) that this function decays at the 
special transition as 

{ r-^d-2+rnie-i)+n,] for r„ = 0, 



as well 


as 


the surface critical indices 








= ri + ril 




and 














P? 


= \v[d- 


2 + ,- 



r 



-[d-2+m(6/-l)+»7||]/e 



for rp = 



in the limit of large separations r (parallel to the surface). This translates into the 
small-momentum behaviour 

^"^^^^ip(^"-)/«, forp, = 0, ^''^ 

of the momentum-dependent local surface susceptibility xii(p) = G'"'^^(p) (its Fourier 
transform), where p ~ {pa,pp) is the wave-vector conjugate to r = (r^, rp). We have 
suppressed the specification 'sp' on ?7|| in these two equations because analogous results 
hold in the case of the ordinary transition, with ?7|| given by rj?,'^'^, its analogue for the 
ordinary transition (Diehl et al 2003a, 2003b). For m — 0, the above dependences on 
r and p for ra — and pa ^ 0, respectively, reduce to the familiar power laws for the 
case of the special transition at a bulk critical point (Diehl and Dietrich 1981, 1983). 
Henceforth we will refer to this surface transition as the 'critical-point (CP) special 
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transition', calling its tti > analogues we are concerned with here 'Lifshitz-point 
(LP) special transitions'. 

Note also that the scaling forms (49) can be exploited in a straightforward fashion 
to generalize the standard scaling laws of the CP case to to > 0. Consider, for example, 
the surface critical exponent 77^ defined through the decay 

G(i'i) [(r, z), r] ^ z-[d-2+rn(6-i)+n,_] (52) 

at the LP.^ Upon introducing the momentum-dependent layer surface susceptibihty 

Xi(p)= / dzG(i'i)(p;z), (53) 

Jo 

we find from equation (52) the small- momentum behaviourxi(p) ^ p^^^"^ or 
r^ p(')i-2)/e^ depending on whether pa = Q ot pp = 0. The result (49) implies that 
the familiar scaling law 

?7± - i(r7 + r?||) (54) 

holds. Furthermore, let 711, 71, j3a and 7^ denote the critical indices that characterize, 
at (0 = c = 0, the thermal singularities of the surface susceptibilities xii = Xii(O) ^^id 
Xi = Xi(0), the surface excess order parameter 

[G("^i)(r,z)-G(°'i)(r,oo)]dz (55) 

/o 

and the surface excess susceptibility Xs = d?Tis/d/i|;i=Oj its derivative with respect to 
a bulk magnetic field ft,, via Xi'"^ ~ 9t^^^j Xi'"^ ^ 9t'^^i '^s ^ 9t ^^^^ xl'"^ ^ 9t'''' j 
respectively. From equations (49) and (52) we see that the scaling laws 

711 = z^(l-?7||) , 71 == i^(2-77_L) , (56) 

and 

A=/3-i^, 7. ==7 + ;^, (57) 

also remain valid for to, > 0. 

Finally, let us briefly consider the surface energy density ([(/)®(r)]^). In analogy 
to the CP case (Dietrich and Diehl 1981), the renormalization of this quantity 
involves further (additive) counter-terms, and the implied RG equation becomes 
inhomogcneous. We refrain from working out the details here, noting that the form 

([^<B(^)]2)sing^^2-a-.-* (.9,^^=0) (58) 

of its thermal singularity at p =: c = can be inferred from equation (49). 

Next, we turn to the e-expansion of the surface critical exponents of the special 
transition. From equation (15) we see that the 0(e) contribution to the fixed-point 
value ?7* vanishes for all to, > 0, so that 

Vf -V = vl = 0(6^) (m > 0) . (59) 

On the other hand, substitution of the fixed-point values (28) and (30) into equation 
(16) yields 

--1-'?:--^[1 + Ao(to)z;(Ao;to)]6 + 0(62). (60) 

z/ n + 8 

^ This exponent may be defined more generally through the decay of the pair correlation function 
G^^'^'{x, x') at criticality as x' moves away from the surface along a direction perpendicular (or at 
least not parallel) to it while x' is fixed, similarly as in the CP case (Diehl and Dietrich 1981, 1983, 
Diehl 1986). 
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Let us see what happens to these resuhs in the hmit m ^ 0. Equation (60) 
reduces indeed for to ^ to the correct 0(e) result for the CP case, just as the 
e-expansions of the bulk critical exponents (Diehl and Shpot 2000, Shpot and Dichl 
2001) and the surface critical exponents of the LP ordinary transitions (Dichl et al 
2003a, 2003b) do. On the other hand, the 0(e) term of the exponent rjfF of the CP 
special transition is known to be nonzero (Bray and Moore 1977, Diehl and Dietrich 
1981, 1983). Consequently the limits to, ^ of the e-expansion (59) of rjf,^ (as well 
as those of related surface exponents such as /3^^) differ from those of their to, = 
counterparts for the CP special transition. Nevertheless, this is no cause for concern. 
The obvious reason for the apparent discrepancy between the ?n — > limit of our 
result for rjl^ and the established one for m = is that the limits A -^ \*^{m) and 
to ^ do not commute. Taking the limit ?Ti — > last means to let to ^ while A 
follows the root A^(to) on its way to infinity. Since rii[u* , X*^^{m)] is of order e^, so is 
its limit ?Ti ^ 0. 

However, in order to compare with the CP case, the relevant problem to consider 
is the following. Given a generic initial value A(l) = A < oo, what asymptotic 
behaviour does the RG flow for to, ^ yield on large length scales? Noting that 
the running variable X{£) —^ oo, we see that the appropriate limit of rji that matters 
is lim)j_^QQ limm_,o ''7i('w*, A). This implies the replacement of the integral ii(A;TO,) 
by its TO ^ limit ii = 1 in our result (15) for 771 (u. A), whereby it reduces to the 
correct one-loop expression for the CP case. Upon insertion of its value a,t u — u* into 
equation (59), the familiar result for the exponent rhF of the CP special transition is 
recovered to order e. Thus, extrapolated to to = 0, the RG flows obtained here yield 
results in conformity with the CP case, despite the mentioned difference between the 
TO ^- limit of equation (59) and the e-expansion of r]l^ for the CP case. 

Let us try to use the results (59) and (60) to estimate the values of the associated 
surface critical exponents of the three-dimensional ANNNI model. Since we know 
these series expansions merely to first order in e, only crude estimates are possible. 
Clearly, it is desirable to exploit as much as possible the improved knowledge about 
the bulk critical exponents gained in Diehl and Shpot (2000), Shpot and Diehl (2001). 
It appears preferable to focus directly on estimates of the differences ?7m^ — tjot Pl^ — (3 
(both of which are proportional to 77*) and the ratio $/z^, rather than extrapolating the 
expansions to 0(e) of surface critical indices such as rfv^, /3^^ or surface susceptibility 
exponents and $. 

From the result (59) we conclude that the differences rf^y — rj and /3l^ — (3 are 
likely to be small. Accepting the best estimate 77 ~ 0.124 of Shpot and Dichl (2001) 
for the uniaxial scalar case in d — 3 dimensions, we conclude that ih? ~ 0.1. One way 

of estimating I3f is to employ the values v ~ 0.746, 9 = v^/v ~ 0.348/0.746 ~ 0.47 
of Shpot and Diehl (2001) in equation (48), along with the approximation rf^F ~ 
r/ ~ 0.126. This gives Pf ~ 0.22. If the slightly bigger value e ~ 0.487 quoted in 
Diehl (2002) (and found by direct extrapolation of the e-expansion of 9) is utilized 
instead, one obtains /3^^ ~ 0.23. The agreement with Plcimling's Monte Carlo estimate 
/3^^ = 0.23(1) is impressive. Another possibility is to start from the second part of 
equation (48) and utilize the estimate (f^ ~ (3. The best estimate of Shpot and Dichl 
(2001) for /3 was /3 ~ 0.246. The most recent Monte Carlo result /3 = 0.238 ± 0.005 
(Pleimling 2002, 2004, Pleimling and Hcnkel 2001) is in conformity with this and lends 
support to the approximation j3 ~ j3^^ . 
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Unfortunately, we are not aware of any estimates of the surface crossover exponent 
^ at d — 3. If we take again m = n = I and boldly set e = 3/2 in the expansion 
(60), truncated at 0(e), we obtain ^/ly ^ | x 1.65816 x | ~ 0.83. Utilizing once 
more the value v ~ 0.746 then yields $ ~ 0.62. We must caution the reader, however, 
to take these estimates cum grano salis: the e-expansion to first order does not in 
general give reliable results. Moreover, in the CP case (vn = 0), the e^ term of the 
crossover exponent $ is known to have a rather large coefficient (Diehl and Dietrich 
1981, 1983), a fact which makes it difficult to obtain precise estimates even from the 
e-expansion to second order. Borel-Pade estimates based on the massive field theory 
approach at fixed d ~ 3 (Diehl and Shpot 1994, 1998) give values ^d=3,m=o,n=i — 0.54 
considerably smaller than the original one (~ 0.68) found in Diehl and Dietrich (1991, 
1983), Diehl (1986) by evaluating the e-expansion to second order at d = 3, though in 
much better agreement with Monte Carlo results (Landau and Binder 1980, Ruge et 
al 1992, 1993). 

In extrapolating our e-expansion results for the scalar uniaxial case m — n — 1 
to three dimensions, we took it for granted that an LP special transition is possible at 
d — i. Clearly, a transition to a surface phase with long-range order should be possible 
at temperatures Tg higher than the line of bulk critical temperatures Tc \i m = n = 1. 
Since Tg can be varied through a change of surface interaction constants, one expects 
that for appropriately fine-tuned surface enhancement Tg can become equal to the LP 
temperature, so that a special LP transition ought to be possible at d = 3, as also 
recent Monte Carlo results (Pleimling 2002) indicate. 

For other values of m and n extrapolations to d = 3 make little sense. Since we 
precluded in the Hamiltonian (l)-(3) any (bulk and surface) terms breaking its 0{n) 
symmetry, a bulk-disordered surface phase with long-range order — and hence a special 
transition of the kind considered above — cannot occur for rf = 3 if n > 1 because it 
would require the breaking of a continuous symmetry in an effectively two-dimensional 
system. 

It is tempting to anticipate the possibility of a special transition for the biaxial 
scalar case m = 2, n = 1. However, in the case of a multi-axial LP, one generically 
expects contributions to the Hamiltonian that break its rotational invariance in the 
(m > l)-dimensional subspace of a-coordinates, such as bulk terms of cubic symmetry, 
X]a=i('5a'^)^i and similar space anisotropics of further reduced symmetry. According 
to a recent two- loop RG analysis (Diehl et al 2003c), such anisotropics are relevant, 
at least for small e > 0. Unfortunately, no stable new (anisotropic) fixed point 
could be found. To our knowledge, Monte Carlo investigations of appropriate biaxial 
generalization of the d = i ANNNI model have not yet been carried out. Thus, 
whenever such anisotropics cannot be ruled out, it is presently unclear whether a 
biaxial bulk LP exists, in particular in three dimensions. The clarification of this issue 
is beyond the scope of the present paper. Note also, that in order to account for 
the absence of rotational symmetry in the m-dimensional subspace, we would have to 
generalize the surface part Ci of the Hamiltonian by allowing its derivative term to 
become anisotropic as well. 

4. RG improved perturbation theory 

In the previous section we exploited the RG equations (10) to deduce the scaling forms 
(49) of the (A^+M)-point cumulants (4). A tacit assumption underlying this derivation 
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is that the renormahzcd cumulants Gron are well-behaved at the fixed point V*. 
More precisely, these functions were assumed to approach for e > finite and nonzero 
limits, as A and u approach their respective fixed-point values A^ and u* . This can 
and should be checked within the framework of RG improved perturbation theory. 
Furthermore, other conditions exist to which this statement applies just as much. 
For example, consistency of the scaling forms for gr y^ with their counterparts for 
Qr = imposes constraints on the asymptotic (/T--dependence of the scaling function for 
Qr ^ 0. Further constraints concern the dependence on gc] they arise from matching 
requirements between the (/c-dependent scaling forms and the asymptotic behaviour 
at those transitions to which the special transition crosses over, namely the ordinary, 
surface or extraordinary transition. 

Our aim here is the explicit verification of some of these properties by means of 
RG improved perturbation theory to one-loop order. For the sake of simplicity, we will 
restrict ourselves to the illustrative case of the (renormalized) surface susceptibility 
XiT" (p)- We begin by considering this quantity at t = p = for general nonzero 
momentum p = (pajPp) G M as function oi gc- Upon introducing the dimensionless 
momenta 

P = (pppp/f^'f , P - i'y'^'p.pJt^Y^' , (61) 

we can write its scaling form as 

x[T\p) « ^,-'p^7-' s{Pp-\9cp-^''') (62) 

= ^-lp(''^-l)/%(pp-^5,p-*/-), (63) 

where the scaling functions $ and S are related to each other via 

*(P,c) = p(^"''r)/''s(P,c) . (64) 

Consistency with the scaling forms that hold if cither p = or P = requires the 
limiting behaviour 

So(c) forP^O, 

'*oc(c) forP 



^(^•='Hp(„;.-.)/» - '«^' 



Next we consider the limits p ^ and P ^ for gc > 0. The limiting dependences 

ord__-| 

on p and P must be in conformity with the leading infrared singularities ~ p''ii 
and ^ P"ii ~ •'' at the ordinary transition. Hence we anticipate that S varies for 

c ^ oo as 

s(p,c) « c-'^(i-''r)/*[xn+yiP'c-2-''/*+...]+c-'^("r'-''r)/*Soo(p)+...,(66) 

where Xn and Yi are constants, while Soc(P) behaves as 



"00,0 



S„,(0) forP^O, 



^-<^'^^ .or-.)/. f„p^eo. ""' 



^00,00 I 



The latter properties ensure that the term of equation (66) following the square 

brackets produces the required powers p**" ~^ and P'-''ii ~^'' . To understand the 
terms of equation (66) involving the square brackets, one should recall that xn remains 
finite at the ordinary transition. Thus a vaovaBnixmi- independent contribution to 
X-l^"]^" (0) must exist, apart from additional ones with analytic momentum dependence. 
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The contribution oc Xn accounts for the former; except for a contribution proportional 
to P^, terms of the latter analytic kind have been suppressed, along with momentum- 
dependent subleading corrections. 

In Appendix B the renormalized function x^'j^" is computed to one-loop order 
for T = p = and general values of c > and A > 0. The result can be written as 

n ' „ u n + 2 

2 

,•(0,1) , o^^ A(i,l f2 „l^ \\ , r^l.\\ , r^/„.2 



(rcn) / \ 

Xii '{p) 



= K„ + c+AP^ + ^ '-^-^ \ [{Kp + c) «i + c (l + A «i)] ln(2Kp) - c\ i^^"^^ 



- {np + c) i\ 



2-rA{P/^,c/p,X)+0ie) 
p ^ ' 



0(u' 



where 



P4 



while the function A is an average of the form (21): 



-4(P,c,A) 



In [c/(l + At)VTTP^ 



(l + At)2[c-(l + At)\/TTP4l 



(68) 



(69) 



(70) 



,•(0,1) 



Further, i"^^'^' and i\ ' denote the partial derivatives dii/de and d^ii/dedX, 



respectively. Just as ii and {[, these are to be evaluated at e = in equation (68). 

From the perturbative result (68), the universal scaling function S(P,c) must 
follow up to non-universal metric factors upon setting the coupling constants u and 
A to their fixed-point values (28) and (30)+. Recalling the 0(e) expressions (59) and 
(60) for ?7y'^ and $/j^, one easily verifies that xj^"]"" takes indeed the scaling form (62), 
where the scaling function is given by 



S(P,c) 



with 



and 



a„(e) VT+P^ + bn{e) c + Ao P" + ^^ e|2 c^ ^(P, c, Ao) 

n + 8 L 

+ c[l + Ao^'l(Ao;m)] In [2^1 + P^] } + 0(6^) 



a„(e) = 1 - e 



bn{e) = 1 



n + 2 


n + 8 


n + 2 



z^^)(Ao,0;m) 



if-'\Xo,0;m) + Xoi']'-'\Xo,0;m) 



(71) 



(72) 



(73) 



Obviously the limiting forms (65) for small and large P hold. Checking the 
consistency with the critical behaviour at the ordinary transition is less trivial 
because the e-expanded version (71) of the scaling function docs not yet contain 
the anticipated asymptotic power laws (65) and (66) in exponentiated form. To 



proceed, we consider the behaviour of ^ as c 



). The leading contributions 
P** of this order cancel. The 



are of order c In c. One finds that the terms ex In \'l 

remaining (momentum-independent) terms ^ cine can be combined with part of 

the contributions ^ c to identify the power c^^'^" to 0(e). The associated exponent 

"*" As non-universal factors we have, in this case, the overaU amplitude of S(P, c) and the two metric 
factors associated with a and gc = c, respectively. 
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satisfies 1 — 77* = v{l — rfJ')/^ + O(e^). Hence tliis contribution gives rise to the first 
term in equation (66). For its amplitude Xn we obtain tlie e-expansion 

11 „ I o 



•("•1)m o ^ I \ •(1-1)^ o ^ I o /^^(-"^ + ^o^)\ 
i\ '(Ao,0;m) + Ao^l '(Ao,0;m) + 2( -rr-— t-tt ) 

-[l + Ao^'l(Ao;™)]ln2|+0(e2). (74) 

Turning to the contributions of order Inc and c", we note that the former can be 
combined with part of the latter to obtain the term e^jif n ln(c/K), where k = Vl + P^, 
and we employed the fact that ((1 + Ao)~^)o,m = («i(Ao;m) + (l)o,m)/2 = 1/2. This 
logarithm can be cast in the form k [(c/k)'^'^"'"'"^^/'^"^*' — 1] + O(e^). Recalling the 
e-expansion 

2-e^ + 0(62), (75) 

" n + 8 

one easily verifies that the term k (c/k)'^*^"+^)/("+®' produces a contribution to S 
consistent with the form of the term in the second line of equation (66), with 

Soc(P) = -(1 + P^)'''"'"''^'' [1 + 0(6)]. (76) 

Furthermore, the term ex Aq P^ in equation (71) yields an analytic contribution ex P^. 
Its form matches the one ex Yi in equation (66), where the coefficient is given by 

Yi = -Ao + 0(e) . (77) 

In summary, we find that our result (71) for the scaling function S is also in 
conformity with the limiting forms (66) and (67). This means that the asymptotic 
momentum singularities one obtains at the ordinary transition from the c-dependent 
analysis here, to first order in e, agree with those found in Diehl et al (2003a) for the 
case c = 00 with Dirichlet boundary conditions. While expected, this manifestation of 
universality is gratifying in that it corroborates our identification of the fixed points 
Kp and Vl,^. 

5. Summary and conclusions 

We have investigated the special surface transition that occurs at an ?Ti-axial bulk LP 
in the case where the m (potential) modulation axes are parallel to the surface. Our 
RG analysis vi\d = A-\-'!-!^— e dimensions corroborates the observation of Diehl et al 
(2003a, 2003b) that a square gradient term of the form shown in equation (3) must 
be included in the boundary part Li of the Hamiltonian, in addition to its usual (jp' 
contribution. 

The RG fixed point associated with the special transition is located at a nontrivial 
value A!j_ = 0{e^) of the associated renormalized surface coupling constant A, zero 
(renormalized) surface enhancement c = (see equation (29)) and the usual fixed- 
point value u* of order e of the bulk interaction constant. As a direct consequence 
of the fact that \\ = 0(e°), the surface crossover exponent $ becomes m dependent 
already at first order in e, whereas this happens for the bulk critical indices and the 
surface critical exponents of the ordinary transition only at second order in e. 

It would be natural to anticipate the same behaviour as that of <& for the surface 
correlation exponent rftF , the second independent surface critical index of the special 
transition. Yet this is not the case because the expansion to order e of r/^^ turns out 
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to vanish, whereas the contribution of order e of its ?7i = analogue is known to be 
nonzero (Diehl and Dietrich 1981, 1983), as is its O(e^) term. Thus the hmit m ^ of 
the e-expansion of rhF differs from the estabhshed m — Q result for the CP case. This 
behaviour is exceptional. For once, the e-expansions of all those hulk critical indices 
that retain their physical meaning for to = turn over into their m = Q counterparts 
for the CP case in the limit m — > 0, as has been explicitly verified to order e^ (Shpot 
and Diehl 2001). Second, the same applies to the surface critical exponents of the 
ordinary transition (Diehl et al 2003a, 2003b). Moreover, even our 0(e) result for $ 
implied by equation (60) is fully consistent with the e-expansion of $m=o (Diehl and 
Dietrich 1981, 1983). 

As we have seen, the origin of the discrepancy between the e-expansions of 
limm_,o r/^tF and its ?Ti = analogue for the CP case can be traced back to the non- 
commutativity of the limits A — > \\{m) and to, ^ of the RG function rii(u* , A). It 
is gratifying that we could recover the correct 0(e) term of rjl^ for the CP case by 
analysing the RG flow directly for to ^ 0. 

The above, not necessarily expected, results are a clear message that further 
surprises might well be encountered when the surface orientation is taken to be 
perpendicular. Since in this case the distance z from the surface scales naively as 
y^-i/2 (rather than as the inverse of the momentum scale n), more surface monomials 
exist whose coupling constants have nonnegative momentum dimensions for e > 0, 
and hence are potentially dangerous. This makes the analysis of multi-critical surface 
transitions such as the special one even more interesting and challenging. 

The series expansions of the surface critical exponents of the LP special transition 
determined here arc restricted to first order in e. Unless their knowledge can be 
combined with information from other sources, one certainly cannot hope to get 
numerically accurate estimates of the surface critical exponents at d ~ 3. The 
additional information we could benefit from was the e-expansions of the bulk critical 
exponents to O(e^) and previous field-theory and Monte Carlo estimates for d — 3 
(Shpot and Diehl 2001, Pleimling and Henkel 2001). Our result (59), which means 
that Pl^ = /3 -(- O(e^), is in accordance with the small difference between the values of 
Pl^ and the bulk exponent /3 obtained in Pleimling (2002) via Monte Carlo simulations 
of the three-dimensional ANNNI model. 

We see no principal obstacle to extending the present analysis to second order in e. 
However, the effort required appears to be considerably greater than was necessary for 
our — already quite demanding — analysis of the ordinary transition (Diehl et al 2003a, 
2003b). Furthermore, the residues of some pole terms one must determine seem to be 
expressible only in terms of multi-dimensional integrals which must be computed by 
numerical means. 

In order to corroborate our identification of the fixed points (29) we also computed 
the scaling functions (71) of the momentum-dependent local surface susceptibility 
Xiiip) ^t the Lifshitz point r = /? = 0, to first order in e, and confirmed that the 
momentum singularities at both the special as well as the ordinary transition have 
the predicted power-law forms, consistent with our own RG results here and those of 
Diehl et al (2003a) for the ordinary transition. 

We close with three remarks. 

(i) There are several worthwhile goals which future Monte Carlo simulations of 
the ANNNI model might accomplish. An obvious one is obtaining an accurate 
estimate of the surface crossover exponent $. A second — equally important and 
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not unrelated — task is a careful check of the predicted scaling behaviour at the 
special transition. We believe that appropriate (non)linear scaling fields g^, gr etc 
[see equations (9)] should be introduced when performing the required analysis 
of the data. 

(ii) While mathematically rigorous results on the surface phase diagram of the three- 
dimensional semi-infinite Ising model can be found in the literature (Frohlich and 
Pfister 1987, Pfister and Penrose 1988), we are not aware of similar work on the 
semi-infinite ANNNI model. Clearly, rigorous results both on bulk and surface 
phase diagrams involving bulk LP could be very valuable. 

(iii) Finally, let us emphasize the need for careful experimental studies of bulk and 
surface critical behaviour at LP. Owing to the enormous advances in experimental 
technology and our better theoretical understanding, tests much more stringent 
than the (more than 20-year-old) experimental study (Shapira et al 1981) of bulk 
critical behaviour at the uniaxial LP of MnP ought to be possible today. The 
situation is even worse for surface critical behaviour at LP, since experimental 
investigations do not yet exist apparently. 
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Appendix A. Calculation of Ac 



In this appendix we derive the one-loop result (13) for the renormalization function 
Ac- To this end, we consider the surface cumulant Gron (p) for t = c = and generic 
A > and p > 0. The free propagator of the rcnormalized theory may be gleaned 
from equation (7) of Dichl et al (2003a) and the rc-parametrizations (5) and (6). We 
express it in terms of the dimcnsionless momenta p and P introduced in equation (61). 
For T = c = 0, it then becomes 



G(p;zi,Z2) = 
in the pz representation, with 



1 



^Lirhjj 



^ — Kp\Zl-Z2\ 



K„ + AP2 



p4 _^ pp2 



pz 



(A.l) 



(A.2) 



Upon taking into account all counter-terms that contribute to first order in u, we 
obtain 



«p + (A + Pa)P2 

/•OO 

dz [G{p; 0, ; 



Q 



- ArfJ, p - Acpp5{z) 



0{u^ 



(A.3) 



where the lower limit 0— of the z integration serves to ensure that J dz5{z) = 1. 
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The tadpole graph displayed in this equation yields two contributions: one 
corresponding to the translation invariant part of G — its 'bulk' part Gb — and a second 
one resulting from its remaining 'surface' part Gg- Since 

Gb(---)=y(^^prT^^==^;;,7^^ + 0(O, (A.4) 

b 

the pole term of (^ (the former) gets cancelled by the subtraction ex A^- in the square 
brackets, for our choice of Ar, made in accordance with equation (17) of Diehl et al 
(2003a). 

The graph ^ yields pole terms of the form 5'{z)/e and 5{z)j€. To see this, note 
that upon making the change of variables p -^ pjl^iz and P -^ P / ^2pz^ the surface 
part Gs{x, x) can be written as 

Gs(aJ, x) = a-"/4(2z)^-Ve(A*^) (A.5) 



with 



where 



^pd^^P K,{p,P,z)-XP^ e-«(P'-P^^) 



d — 'm—l' 

' (27r)'^-i k{p,P,z) + XP^ 2Kip,P,z) 



k{p, P, z) = y p2 + P4 + 2zpP^ . (A.7) 

The pole terms of the generalized function (A.5) originate from the contributions that 
behave as z'^"^ and z'^~^ near z = 0. Utilizing the well-known Laurent expansion 
z'^~'' = {—l)''S^''^z)/k\ + 0(e) of such generalized functions (Gelf'and and Shilov 
1964) gives 

(2i)-2 ;^(^) = --IMs'iz) + IM. 5(2) + 0(e) . (A.8) 

The Taylor coefficients /(O) and /'(O) can both be expressed in terms of the function 
Zi and its derivative i[. A straightforward calculation yields 

/o(0) = 2F„^on(A;m) (A.9) 

and 

'l-ii{X;m)- Xi[{X;m) ^^^[(to + 2)/4)]" 



/^(0)--4pF„ 



(A.IO) 



2A 2r(m/4) 

s 

The uv singularity of the graph (^ implied by the pole term ex 6{z) in equation (A.8) 
is momentum independent; it must be cancelled by a renormalization of the surface 
enhancement c. Since we set the renormalized quantity c to zero, this condition 
immediately gives the renormalization function Ac to one-loop order. The result (13) 
can be read off from the above equations. The pole term ex S'{z) in equation (A.8) 
looks like a counter-term of the form J^ (pdn<p, which we did not introduce, knowing 
that it can be transformed by means of the boundary condition 

9„0==(c-Aa„a„)0 (A.ii) 

into surface counter-terms we included, such as the one ex jm{da4>)^ (Diehl et al 
2003a). In fact, using the one- loop results 

Zi(u, A, e) - 1 = -X-'P,{u, A, e) + 0{u') = ^ !i(^!^ + o(z.2) (A.12) 

of this reference, one easily convinces oneself that the pole terms produced by the 
renormalization functions Zi and P\ of the first term in equation (A. 3) cancel those 
originating from the singularity ex S'{z) of Gs{x, x). 
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Appendix B. One-loop calculation of Xiiip) 

In this appendix we outline the computation of the renorniahzed surface susceptibihty 
Xn " (p) to one-loop order for t = p = and general values c > and A > 0. The 

b 

contribution of the bulk part (^ to the one-loop term vanishes for t = p = (in 
dimensional regularization), just as the terms proportional to At and Ac in equation 
(A. 3) do. In the remaining one-loop integral we perform both the integration over z 
as well as the angular integrations, obtaining 



" ^Z, 



(rcn) / ■ 

Xii (p 



«p + (A + Pa)P^ + ^cC +^^ ^'^7'^^" J + 0(u^) (B.l) 



oo 



with 

/•ex 

J= dpip^-"'-'-' I dPiA""' '— — '-^^ ^- (B.2) 

Jo Jo 4Kpj(Kp + fi;pJ(Kpj + c + APf) 

Here p and P denote again the dimensionless momenta (61), while Kp, defined by 
equation (69), is the p — analogue of the quantity (A. 2). The renormalization factor 
Zc may be gleaned from equation (22) of (Dichl et al 2003a): expressed in terms of 
i'l, it reads 

Z,(u,A,e) = l + ^-^ ^^ ' '^ +0{u^). (B.3) 

Finally, Zi and Pa are given by equation (A. 12). 

In the double integral J we introduce new variables s and t via P = yt/s and 
p = s^^^/l — t^, whereby J becomes 

j='-[\tt^--y^[i-t^f—^^^'' rds/'^'^^y--^") 

8io ^ ' Jo (1 + sKp) (1 + cs + tA) 

^^{e^/ r^ s-'+' jl -cs-tX) \ 

8 \io ■ (l + s«p)(l + cs + tA)/^,„ ' ^ • ^ 

where {■)e,m and A^e,™ are the average and the normalization factor introduced 
by equations (21) and (22), respectively. The integral /„ ds can be performed 
analytically and its Laurent series coefficients of order e~^ and e'' determined. One 
thus arrives at 

uK,.rr.-iKmj^ -^|[6-i-ln(2«p)] [«p^l(A,e;™) + 2c((l + At)-2)^ 



2F 2 



ln[c/Kp(l + At)] 



\(l + At)2[(l + At)Kp-c]/o„^ J 

We now substitute this result together with the above expressions for the 
renormalization functions Zi, Zc and P\ into equation (B.l) and use the fact that 
the average in the first line of equation (B.5) can be written as 

((1 + At)-2)^^^ = i [1 + zi(A,e;m) + Xz[{X,e;m)] . (B.6) 

The poles in e are found to cancel, and one arrives at the result (68). 
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